We propose a model for rotating current patterns 
within radiopulsar polar cap accelerators which has obser- 
vational consequences that mimic those which have been 
attributed to neutron star precession. The model is a sim- 
ple extension of a commonly used one for the origin of 
the "drifting subpulses" often observed within the pulse 
envelope of radiopulsars. The new model's current pat- 
tern rotation period (with respect to the neutron star) is 
estimated to be of order a year. Associated with that ro- 
tation are small oscillations in spin-down torque, pulse ar- 
rival time, and radiobeam direction with this same period. 
These have estimated magnitudes which support a pos- 
sible reinterpretation of precession "observations" which 
- - , could resolve the problem of obtaining the required pre- 
cession parameters with canonical neutron star models. 

o 

{S| 1 Key words, pulsars: general - stars: neutron - pulsars 

o 
O 



> 

00 

m 
O 

o 

S3 

Pl,: 

6 



X 



Astronomy & Astrophysics manuscript no. ruderman 


February 5, 2008 


(DOI: will be inserted by hand later) 





Mimicking neutron star precession by polar cap current-pattern 

drifting 

M. Ruderman 1,2 and J. Gil 2 ' 3 

1 Department of Physics and Columbia Astrophysics Laboratory, Columbia University, New York, NY 
e-mail: mar@astro.columbia.edu 

2 Institute of Astronomy, Cambridge, U.K. 

3 Institute of Astronomy, University of Zielona Cora, Lubuska 2, 65-265, Zielona Gora, Poland 
e-mail: jag@astro.ia.uz.zgora.pl 



Abstract. 



1. Introduction 

For a considerable number of pulsars the residuals (At) 
between observed pulse arrival times and a smooth expec- 
tation based upon average period P, P, and P can be fit 
to a single oscillation 

At = aSin (i?) (1) 

and its overtones (cf Table 1). The physical reality and 
significance of this residual oscillation is particularly com- 
pelling for PSR 1828-11 (Stairs, Lyne, & Shemar IMTH)! 
whose observed time variation in pulse shape is very 
strongly correlated with an oscillating component in P(t), 
with the same period P p . This correlated pulse-shape com- 
ponent has discouraged an explanation of the At oscilla- 
tion as neutron star (NS) movement in response to that of 
planetary companions. Instead it has generally been con- 
sidered as strong evidence that this neutron star has a 
sustained precession with period P p of about 3 years. 

A long-lasting precession period of such magnitude is 
compatible with pulsar models which describe a neutron 

Table 1. Period (P p ) and amplitudes (a) for the oscillat- 
ing component in the time residuals between observations 
and expectations for several radiopulsars. 



Pulsar 


P p (days) 


a{10^ sec) 


Ref 


Crab 


10 a 


1 


LPS88 


B1828-11 


10 3 


2 


SLSOO 


B1642-03 


2 x 10 3 


2 


SLU01 


B1929+10 


1 x 10 3 


1 


U03 


B1557-50 


2 x 10 3 


0.5 


CUO03 


RX J0720 


3 x 10 3 


10"(?) 


HTD06 



star as a rigid body (Agkiin, Link and Wasserman 2006; 
this paper also contains an excellent well-referenced survey 
of previous work on NS precession.) However, the preces- 
sion parameters inferred from observations are much too 
long to be compatible with the canonical models of a NS 
core containing a neutron-superfluid and a proton- Typell- 
superconductor threaded by a strong magnetic field link- 
ing the core, crust, and external magnetospherc (Link 
2003; see also Shaham 1977, Sedrakian, Wasserman and 
Cordes 1999, and a different view by Alpar 2005). In 
a proton- Typell-superconductor, magnetic field is orga- 
nized into quantized flux-tubes which interact strongly 
with coexisting quantized vortex-lines of the spinning neu- 
tron supernuid. If, in a less standard model, the protons 
were to form a Type I-superconductor, a "mixed state" 
would appear in which the magnetic field becomes strong 
enough in some regions to quench superconductivity there 
but vanishes elsewhere. According to Sedrakian (2005) the 
drag on vortex-lines moving inside the very different re- 
gions of such a "mixed state" can be small enough to 
be compatible with long-period NS precession. (However, 
special consideration of expected pinning or pile- up of vor- 
tices at very numerous interfaces between magnetized and 
unmagnetized regions still seems needed). 

The canonical model for a NS, with its strong interac- 
tion between a core's moving superfluid neutron vortex- 
lines and the magnetic field out through the NS's con- 
ducting crust, agrees well with a number of different 
kinds of pulsar observations: dipole magnetic field evo- 
lution in spinning-down pulsars, magnetic field structures 
of millisecond pulsars, Crab-like "glitches" , giant Vela- like 
"glitches" (Ruderman 2005, 2006). It seems problematic 
whether a non-canonical model with a very much weaker 
interaction between core vortex-lines and surface magnetic 
should have similar consequences. 
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We propose below models which avoid real and po- 
tential conflicts between canonical models of neutron star 
interiors and the Table 1 " precession parameters" . In them 
current-patterns in the pulsar polar caps (and magneto- 
sphere) rotate in a polar cap -B-field which remains fixed 
in a non-precessing neutron star. The estimated current- 
pattern rotation period (P p ) and the amplitude (a) of the 
pulse time of arrival oscillation caused by the current ro- 
tation's associated oscillation in NS spin-down torque are 
both similar to those in Tabled 

In Section 2 we first review an old model (Ruderman 
& Sutherland I1975f) for the common rapidly "drifting 
subpulses" observed in many radiopulsars (Backer 119731 
Weltevrede, Edwards, & Stappers (2005 ) . In it these are a 
consequence of E x Be/ B 2 drifting of patterns of "sparks" 
in polar cap accelerators. The estimated and observed cir- 
culation times of the rotation pattern of these sparks in the 
polar cap is tens of seconds. We then consider, in Section 3, 
models for a family of related "drifting spark" models 
which give the very different magnitudes of Table ^ 

2. A rotating "spark" model for canonical drifting 
subpulses 

Deshpande & Rankin 19991 l2()()"Hl have made the most de- 
tailed observations and analysis of narrow subpulses mov- 
ing through pulse envelopes in radiopulsars. Particularly 
striking is their deduction of the arrangement and mo- 
tion relative to the NS of very many narrow radioemis- 
sion beams in a single pulsar. In their analysis of PSR 
0943+10 they deduced 20 separate sources fixed on a 
moving "carousel wheel" rotating in the polar cap's mag- 
netic field with a period ~ 40 sec. Such a moving pat- 
tern has been proposed to come from a similar struc- 
ture in the current flow along the open field line bun- 
dle of magnetic field lines connecting a pulsar's polar cap 
to its light cylinder. Electric fields are needed to main- 
tain such outflow. Without such fields, in the corotating 
frame the inner magnetosphere's electric field (E) and ve- 
locity (v) could be everywhere zero. In the laboratory 
frame v = O x r, E = v x B/c, and the inner magne- 
tosphere would have the Goldreich- Julian charge density 
Pgj — V • E/47T = 1~2 • B/27rc at all points. Because of cur- 
rent flow through the light cylinder (and especially when 
\Sl ■ B| increases with altitude above the polar cap), more 
charge must continually be pulled from the polar cap to 
maintain that flow than can be supplied easily. The polar 
cap surface may not be an adequate source for the needed 
particles because of ion inertia, or perhaps, in some cases, 
because surface ions are bound at a condensed matter in- 
terface. This extreme case is sketched in Fig. 1. (If ions 
can be pulled out but their initial velocities are suppressed 
because of ion inertia the positive excess (space charge) 
would be inside the accelerator, near its lower outerbound- 
ary. The model's consequences would not be qualitatively 
changed. If charge from the NS surface is unavailable, the 
p = gap will grow until the additional potential across it 
extends to a height h and e± pair production can be sus- 
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Fig. 1. Sketch of an idealized polar cap accelerator. The 
magnetosphere plasma on all closed magnetic field lines 
corotates with the neutron star (hatched regions) sepa- 
rated by an empty "gap" from a detached inner magne- 
tosphere plasma (DIP). The positive charge (+) which 
would have existed in the gap region before the acceler- 
ator formed (Goldreich & Julian charge density pgj) is 
now kept from reentering the gap by the stellar surface 
(e.g. by ion binding to a condensed matter surface) and by 
the surrounding open-closed magnetic field line boundary 
where E • B = 0. (Ion inertia alone would have a quali- 
tatively similar consequence with ion pile-up on the gap 
side of the surface boundary). These form an outer sur- 
face whose potential is unchanged by the insertion of an 
accelerator (V = 0). The inner boundary of the accelera- 
tor has a constant V (~ 10 12 volts) where it is just able 
to sustain the e± production needed to limit accelerator 
growth. Effectively, the accelerator contains an additional 
negative charge density p = —pej- The inner boundary 
shape is just that required to make E • B = at the inner 
boundary of the accelerator so there is no charged layer 
there and no change of E within the DIP because of the 
accelerator. In this model the DIP corotates with the neu- 
tron star. The minus signs (-) indicate the absence of the 
positive Goldreich & Julian charge density need to keep 
V = everywhere. 



tained. These pairs are produced by very energetic curva- 
ture 7-rays from accelerated e + and e _ in the accelerator. 
Copious pair production continues far above the acceler- 
ator from those e + (e~) which flow out from it into that 
region. Then because of the consequent high conductiv- 
ity there, E ■ B — in the plasma above the accelerator 
where nearly this same V is maintained from the top of 
the accelerator to near the light cylinder at r ~ r; c = c/Q. 
The indicated potentials and charges in Fig. 1 refer only 
to the differences from those where p — pej everywhere. 
The NS and all plasma on closed field lines (where p still 
equals pgj) still corotate. They are bounded by the V = 
surface. The upper surface of the empty accelerator zone 
has a fixed potential V and a position and shape such 
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that E • B = as well as E tangential vanishing on it in 
the corotating frame. Then, except at the corners where 
the accelerator's outer boundary bends suddenly upwards 
along the last open field line, the accelerator thickness 



( v 



\2irp, 



GJ 



1/2 



(2) 



The entire detached inner magnetosphere plasma 
(DIP) above the accelerator is copiously filled with e 
pairs and it is assumed to have this same V throughout 
it. It does not have any change in its initial corotation 
(E x B = 0) with the NS. In this model the corotation is 
broken only within the accelerator where a charged par- 
ticle will move with respect to the NS with an angular 
speed about a polar cap axis 



E x B 



B 2 



c 



Vc 
hBr\ 



IE x B| 



(3) 



|E x B|, the sin of the angle between accelerator E and 
NS B is at the model polar caps axis (r± = 0) and ~ 1 



at the corner where 



(radius of polar cap) and in 



the slot above it (see Fig. 1 for illustration). 

With the very crude extrapolation \E X B\ ~ r±/r pc , 
equation [3] becomes 



Vc 2 i 



hB d QR 3 ' 



(4) 



This estimate for tt p is not sensitive to spark location 
("slot", near the periphery r± ~ r pc or much closer to 
the r± = axis), and we indicate locations of 4 exem- 
plary sparks in Fig. 1. If the current flow pattern through 
the accelerator consists of particle flow in both directions 
bunched into "sparks" separated by a distance ~ h, a sta- 
ble pattern of such sparks (see Beskin 1982, Gil & Sendyk 
2000) could rotate at the equation J3J rate (e.g. returning 
spark particles heat a surface spot which determines the 
next preferred emission spot, etc.). 

For the pulsar parameters of B1828— 11, P = 0.4 s, 
dipole Bd = 5 x 10 12 G, "characteristic spin-down age" 
10 5 yrs, polar cap radius 2 x 10 4 cm and accelerator thick- 
ness ft~6x 10 3 cm (chosen quite arbitrarily to give some 
partial filling of the otherwise empty accelerator gap be- 
cause of ion flow into it; e.g. Gil, Melikidze & Geppert 
2003) and V ~ 10 12 volts, P p = 2tt/Vl p ~ 40 sec and the 
number of circulating sparks on each circulating carousel 
wheel of the sparks pattern 



27rr 



pc 



20. 



(5) 



The good agreement with £l p and N s deduced by 
Deshpande & Rankin 1X3531 for PSR B0943+10 is com- 
forting but far from compelling and very different kinds 
of models have also been proposed for such drifting sub- 
pulse observations in which outer-magnetosphere plasma 
plays an important role (cf. Wright 2003). Weltevrede and 
colleagues (2005) have shown the lack of correlation be- 
tween observed circulation periods of drifting subpulses 



and various other pulsar observables. Moreover, there is 
often great uncertainty about inferring needed polar cap 
accelerator parameters from those observables. [For ex- 
ample, in several pulsars, and probably in a large fraction 
of them, the polar cap magnetic field seems to be very 
much larger (and the polar cap area very much smaller) 
than the surface dipole field inferred from the pulsar's 
spin-down torque]. A relatively robust prediction of the 
drifting polar cap model of equations (4) and (5) comes 
from combining them to give a predicted angular rotation 
speed of the polar cap's spark-determined current pattern 
independent of the choice of polar cap radius or acceler- 
ator height, but depending only on the number of sparks 
in the carousel. For PSR 0943+10, where that number 
has been determined, the agreement between model and 
observed carousel spin-rate is good. In summary, models 
based upon carousels of polar cap sparks, predict carousel 
rotation periods (of order tens of seconds) and numbers of 
sparks which are not inconsistent with drifting subpulse 
observations. 



3. Very long period circulations of polar cap 
currents and their consequences 

The carousels discussed above rotate because accelerated 
particles from a polar cap surface make pairs, one member 
of which returns to the surface at a point displaced from 
the release point of its parent. The next generation repeats 
this process, but begins from the new origin. Carousel ro- 
tation periods of tens of seconds occur because the time 
it takes an electron (positron) from a pair created within 
the accelerator to return to the polar cap surface is the 
same as the time it spends "drifting" in the strong elec- 
tric field of the accelerator. Rotation periods up to 10 6 
times longer than those considered above would result if 
the pairs were first separated far above the polar cap re- 
gion ( for example, in the outer-magnetosphere near the 
light cylinder). The time of flight back to the polar cap 
surface for a backflowing member of that pair would typ- 
ically be huge relative the time during which it feels the 
strong drift fields in the relatively small polar cap accel- 
erator. 

The processes which make pairs within a polar cap ac- 
celerator continue to make them well above the accelera- 
tor. There is a large pair flux flowing up the open field line 
bundle from the DIP out to the light cylinder. Associated 
with it would be a small potential drop, V' of several 
mc 2 /e ~ 10 6 volts, since an electric field is needed continu- 
ally to adjust local net p to pej because of growing f2 ■ B 
along the pair flow. One way of adjusting the net charge 
density (generally a very small fraction of the total pair 
density) is by the reversal in flow direction of relatively 
few electrons (positrons) which then flow back toward the 
NS surface. The needed readjustment in charge density 
could well be accomplished in other ways. In the absence 
of enough quantitative knowledge of the open field line 
electric field along B out to very near the light-cylinder we 
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simply assume some backflow. If this backfiow dominates 
that directly onto the polar cap surface from the e^ 1 pairs 
created and separated within the accelerator, the RHS of 
equations J3} and JSJ) should be decreased by about /iJl/c. 
Then 



- - ; 



nv 

Br r , r 



V 
~BR 



It) 



1/2 



(6) 



To it should be added an additional rotation from V with 



1 ,:- 



V'c 
R 3 Bfl ' 



(7) 



With PSR B1828-11 parameters, (l' p ~ lO" 1 ^ for pe- 
ripheral and interior spark patterns. For slot gap sparks 
the gap accelerators keep about the same Q p all the way 
up the slot so that equations (J2J and JSJ) remain valid for 
them. For PSR B1828-11 equation © gives P p ~ 0.5 yr. 
For the Crab pulsar P p ~ 0.1 yr. 

The current distribution pattern on the polar cap and 
the j x B forces will determine the large Goldreich-Julian 
torque on a pulsar. We assume that associated with a polar 
cap current-pattern rotation should be a small oscillatory 
component with the same period in the pulsar's spin-down 
torque; SQ =oc rJsin(f2 p t). This would give an amplitude 
in equation Q of 



a 
2t s d 



(8) 



with t s d the characteristic spin-down age (P/2P). We note 
how much larger a becomes as P p increases from the 40 s 
of equation (@J to ~ 10 7 sec in equation |JB}. 

A rough estimate of the dimensionless a is the prod- 
uct of the relative area of any expected rotating irregu- 
larity in the rotating current pattern (h 2 /r pc ) times the 
fractional effect on the Goldreich-Julian spin-down torque 
from moving that current from one side of the polar cap 
to the opposite (r pc /R), times the fraction /(~ 1/2)) of 
the spin-down torque which comes from Goldreich-Julian 
polar cap currents. Then a ~ h 2 /2r p R and, for PSR 
B1828-11 a ~ 3 x 10~ 2 sec. For the Crab pulsar with 
the same V, a ~ 0.4 x 10~ 2 sec. The At of equation JTJ 
oc 2? 3 f2 -1 / 2 so that predicted changes in At are very much 
less than the differences of 10 2 between the Crab and PSR 
B1828's tgd 1 - An exceptionally large amplitude would be 
expected and is observed for RXJ0720. 



1 With realistic polar cap B 7-rays produced near a NS by a 
polar cap accelerator will be gravitationally bent to become a 
strong source of e* creation in large parts of the nearby mag- 
netosphere. Magnetic field structure can have a similar conse- 
quence. These would give e flows on open field lines which are 
not necessarily separated from the polar cap by the accelerator 
potential V. Indeed, such pair flow might even quench parts 
of that accelerator. However the argument for a V' would still 
hold. Q' p would still be present even in the absence of any il p 
for the involved e"/e+. Then for PSR B1828-11 Pp = 4 yrs 
and for the Crab P' = 50 yrs. 



4. Summary 

Our proposed model is a very simple extension of one com- 
monly used to describe the origin of canonical "drifting 
subpulses": E x B drift of localized e± sparks within a 
polar-cap accelerator. It differs only in the time e_(e + ) 
coming back to the stellar surface has spent inside the 
accelerator (h/c, about 10 _7 s) relative to the time spent 
above it (0 — Q , about — 10 _1 s). If the latter time 
= we have the old drifting subpulse model which gives 
carousel circulation times of some 10's of seconds; if it is 
near the maximum time of flight back from near the light 
cylinder, we have the model of section 3. (The choice be- 
tween them may depend on the detailed structure of the 
polar cap's surface magnetic field.) The very rough esti- 
mates of section 4 for P p and At are similar to those of 
Table 1. While all of this is hardly a compelling argu- 
ment for the kind of model proposed here it does seem 
to support further consideration as a not implausible way 
of resolving the problem of how a canonical neutron star 
could sustain the small precession, which has been widely 
suggested as the cause of very long period oscillations in 
timing observables of some pulsars. 
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